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When a cerocene molecule is chemisorbed on metallic substrate, or when an asymmetric double 
dot is hybridized with itinerant electrons, its singlet ground state crosses its lowly excited triplet 
state, leading to a competition between the Zhang-Rice mechanism of singlet-triplet splitting in a 
confined cluster and the Kondo effect (which accompanies the tunneling through quantum dot under 
a Coulomb blockade restriction) . The rich physics of an underscreened 5 = 1 Kondo impurity in the 
presence of low-lying triplet/singlet excitations is exposed. Estimates of the magnetic susceptibility 
and the electric conductance are presented. 



1. The problem of tunneling through a sequence of 
resonance levels was formulated about three decades ago 
[01 . At that stage, little attention was paid to the spe- 
cific structure of the tunnel barrier. Nowadays, novel ex- 
perimental techniques enable the fabrication of artificial 
objects which carry some of the salient features of com- 
plex quantum systems existing in Nature, and to include 
them as resonance barriers in electron tunneling devices. 
They manifest numerous unusual properties and might 
be regarded as important ingredients in future microelec- 
tronics 1^ . Examples are double quantum dot structures 
[Q, atomic and molecular wires and bridges ^j^- 
At the same time, new methods of tunneling mi- 
croscopy make it possible to elucidate the properties of 
single atoms and molecules adsorbed on a surface. The 
combination "nanotip - atom/molecule - substrate" is 
then a quantum system with exceptional resonance fea- 
tures and potential applications 

In the present work we expose the physics of tunnel- 
ing through real and artificial molecules in which there 
is presumably a singlet ground state with an even num- 
ber of electrons, which are spatially separated into two 
groups with different degree of localization. Electrons in 
the first group are responsible for strong correlation ef- 
fects (Coulomb blockade), whereas those in the second 
group are coupled to a metallic reservoir. Hybridization 
with itinerant electrons result in transformation of the 
nonmagnetic (singlet) ground state into a magnetic one. 

Possible real molecules are lanthanocene molecules 
Ln(C8H8)2 with the ions Ln=Ce, Yb in a cage formed by 
TT-bonded carbon atoms 0. In these molecules the elec- 
trons in a strongly correlated /-shell are coupled with 
loosely bound tt electrons. In an analogy with Zhang- 
Rice (ZR) singlet in Cu-0 planes of high-Tc perovskites 
, the ground state of this molecule is a spin singlet 
combination ^Aig{fTr'^) of an f-electron and 7r-orbitals, 
and the energy of the first excited triplet state '^£'2^ is 
rather small (~ 0.5 eV). In the ytterbocene (hole coun- 
terpart of cerocene) the ground state with one /-hole is a 
triplet, and the gap for a singlet excitation is tiny, ~ 0.1 
eV. The fuUerene-like molecules doped with Ce or Yb 
form another family with apparently similar properties. 
In all these systems there is no direct overlap between the 



strongly correlated /-electrons and the metallic reservoir. 
However, these electrons can influence the tunnel proper- 
ties of the molecule via covalent bonding with the outer 
TT-electrons which, in turn, are coupled to the metallic 
reservoir. 

Artificial candidates are double-dot structures (say Di 
and D2) in tunneling contact with each other, but only 
Di is coupled with the metallic leads. The respective gate 
voltages are such that Vgi < Vg2- Coulomb blockade then 
prevents double charging of D2, so it can play the same 
role as 4/ atom in molecular complexes described above. 
The dot Di donates the loosely bound electrons which 
contribute to the tunnel current 

The pertinent physics to be exposed below is that of a 
competition between the ZR mechanism of singlet-triplet 
splitting in a confined cluster and Kondo effect which 
accompanies the tunneling through quantum dot under 
Coulomb blockade confinement jlO|. Usually, tunneling 
through quantum dot containing an even number of elec- 
trons does not display a Kondo resonance due to its spin 
singlet ground state. Analysis of conditions under which 
the singlet ground state changes into a partially screened 
spin-one Kondo state due to hybridization with metallic 
leads is one of the goals of this study. 
2. A simple model which describes this type of molecules 
was considered in |ll], hereafter referred to as a "Fulde 
molecule" (FM). It contains two electrons occupying a 
potential well which is formed by deep and shallow val- 
leys. The Hamiltonian of an isolated FM is 

Here dj^^ creates a dot electron with spin a at valley i ~ 
/, I and spin cr, while the coupling constant V = {di^\df) 
is the inter-well tunneling integral. There are two elec- 
trons in a neutral ground state, and Hcorr = Qnfirif — 
l)/2 is the interaction term responsible for the Coulomb 

blockade of charged states (here Uf = '^/o-'^/o-)- The 
energy difference A — Ei — E f is postulated to exceed the 
overlap integral, /3 = <C 1. Two-electron states |A) 
of the FM are classified as a ground state singlet | S) , low- 
lying triplet exciton |T0), \T±.) and high-energy singlet 



1 



charge-transfer exciton \L). To order they are, 
\S) Ka^\s) -V2l3\ex), 

I™) -^T.d]A-M^ \T±) = rfL4±io>^ 

|i) «a2]|ea;) + V2/3|s), (2) 

where |s) = 75 crc?i-^4^|0), ex) = d\^d\^\Q), and 
= 1 — /9^. In this order, the energy levels i?A are [ pl] |: 

= + £/ - 2^/3, St = e; + £/, = 2(e, + V(3). (3) 

The spin and charge branches of excitation spectrum of 
FM are characterized by rather different energy scales 
Et—Es = S and E^—Es ~ A, respectively. An interplay 
between Kondo triplet excitations (with some character- 
istic energy A^) and ZR triplet excitations is expected 
when S ^ Ak in the regime of Kondo resonance induced 
by tunneling to metallic reservoir fl^ . 

The tunneling problem is encoded in the Anderson 
Hamiltonian which incorporates Hd, together with the 
band Hamiltonian Hf, = X^fccr ^^cla'^f^'^ ^'-'^ elec- 
trons in the leads, and the tunneling term Ht = 
J^ika^i'^kcr'^i'^- Here Cka are operators for lead elec- 
trons and Wi=ij are tunneling matrix elements (as- 
sumed to be k independent but strongly dependent on 
the dot valley quantum number i. It is henceforth 
assumed that Wf = 0). It is convenient to express 
the dot operators di^ in terms of Hubbard operators, 
X^^ = |A)(A|. Here A = S,T,L stands for the neu- 
tral two-electron states ^ , and the index A = la, 3ct 
is reserved for the charged one and three electron states: 
\la) « a\fa)+p\la), |3a) « dy\ex) - ^d] j] ^dl\Q) . 
The tunnel matrix elements in the Hubbard representa- 
tion are given as W^^ = (fctr, A| W^|A), where W is the 
operator responsible for tunneling. The Anderson Hamil- 
tonian then reads: 



Wl^^W, W^'' = ^W, W^^±^W, (5) 

where W = aWi. The energy costs of these transitions 
are 
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H = Y^E^X^^ + Y^e,cl 



Cka 



k<7 



AX 



(4) 



Using the Wigner-Eckart theorem, one can write 
C^\^A, where C;^,, are Clebsh-Gordan coeffi- 
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cients and Ax is the reduced matrix element. In a given 
vector-coupling scheme the tunneling results in the fol- 
lowing transitions: |S),|TO) ^ \l(T,pa); \S),\TO) ^ 
\3a,ka); \T±) ^ |l±,p±); \T±) ^ |3±,fcT). Here pa 
and ka are, respectively, the states with an excess elec- 
tron (and hole) above (below) the Fermi level of the lead. 
Let us focus on the case where the Coulomb blockade 
eliminates the three electron states |3ct) and consider the 
tunnel coupling involving only the states |lcr). The non- 
zero tunnel matrix elements are 



Eip^s = ep - ei + (3V, Eip^r = fp - (6) 
E^ik^s ^ei^AI3V + Q~ e^, E^ik^T = ei + 2f3V + Q - e^, 

where Q « QiV^/iQ - A)^]. 3. We study the interplay 
between the singlet and triplet levels of the double quan- 
tum dot by the renormalization group method following 
the general line of " poor man's scaling" approach to the 
Anderson model The renormalized levels E\ are 

determined by the equations 



dEA/dlnD = Fa/tt. 



(7) 



Here Fa are the tunnel coupling constants, 

Ft = F = ^PoM^, Vs^o?^t, Pg--D-\ (8) 

Integrating (|^) under the conditions E\{Do) — E\, — 
const, we find the scaling invariants E^ which determine 
the scaling trajectories 
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The level £/ is taken to be close to the bottom of the 
conduction band [0, so that scaling does not signifi- 
cantly affect it. It is then subtracted from the energies 
Et and Es- Now we see that the energies Ea decrease 
together with D. Since Ft > F5, the phase trajectory 
Et{D,Tt) should cross that of Es{D,Ts) at a certain 
point. Thus, quite remarkably, there is a crossover from 
singlet to triplet ground state of the FM due to tunnel 
contact with metallic leads. The crossing point can be 
estimated from eqs. (^): nS* « (Ft - Fg) ln(7rD/F), or, 
referring to the bare parameters, the value D of renor- 
malized bandwidth corresponding to this crossing point is 
D = D exp (— ttA/F) . Another important crossing point 
is the energy D = aei (a > 1) where the scaling of e/ also 
stops, the charge fluctuations become irrelevant, and one 
reaches the Schrieffer- Wolff limit where only spin fluctua- 
tions are responsible for scaling of the model Hamiltonian 
p3| . This energy is determined by the equation 



^= (F/7r)exp (^(|e-/-q|)/r). 



(10) 



li D > D, the Schrieffer- Wolff regime is reached after 
crossover from singlet to triplet ground state of the FM, 
and a Kondo type resonant tunneling is feasible. In the 
opposite case there is a singlet ground state and a soft 
triplet exciton (see figure ph. 
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FIG. 1. Scaling trajectories (Eq. |^ demonstrating 
cross-overs from singlet to triplet ground state: Et{D) 
(curve 1), Es{D) for S/F = 0.03,0.1 (curves 2,3, resp.) at 
Ts/Tt = 0.8, Do/T = 10. Inset: S - T transition as a func- 
tion of the level depth ef — at fixed A 

And yet, the 5 = 1/2 Kondo regime is still accessible 
once a properly tuned external magnetic field is applied 
[ p5[ . The novel feature here is that S ^ T crossover can 
be induced by an upward shift of the dot level e; relative 
io ep by a suitable gate voltage (Fig. 1, inset). 

We focus on the physically richer case of triplet solu- 
tion 5 < 0. For a two-electron FM the Schrieffer- Wolff 
transformation |jl^ projects out the states |A, fccr) and 
maps the Hamiltonian H onto an effective Hamiltonian 
H acting in a subspace of two-electron configurations |A) 
and reduced conduction band, H = + + iJ'^-^, 



H 



ST 



tST , 



J'^^ (P • S) . 
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(/i — TO, r±). The local electron operators are defined 

as usual C„ = E/c Cfe^r; S = 2-^/2 J^kk' Haa' 4:a^Ck'a'; f 

are the Pauli matrices. The singlet and triplet states are 
now intermixed, and the spin properties of FM are char- 
acterized by the vector operators S and P in accordance 
with the dynamical symmetry of spin rotator: 

S+ = V2 {X+" + X"-) ,S- = V2 + X-°^ 

S'=X++-X—, P, = -(^^s ' ^S"^ 

P+ ^^/2{X+'^ - X'^-) , P- 



(12) 



[X''^ + X' 

- = V2{X^+ -X-"^), 

These operators obey the moment algebra (i — x, y, z) : 

[P\P^] = ie,,kS\ [P\ S^] = ie^,kP\ S • P = 0. (13) 

and the Casimir operator is + P^ = 3. Surprisingly, 
this special representation of 0(4) played an important 



role in particle physics many years ago 1 14|. The effective 
exchange integrals are 



J^^ = 



^/2A 



(14) 



in which ek is replaced by ej? in the denominators, that 
is, Aa = e_F — €a{D), and eK{D) are the positions of 
the scaled level e; on the flow diagram of Fig. I. Thus, 
the pertinent physics is that of an underscreened Kondo 
impurity ||l6| in the presence of potential scattering and 
low-lying triplet/singlet excitations. A similar model was 
considered recently in Ref. |l^ studying the physics of 
tunneling through a vertical quantum dot in magnetic 
field [|l8|. In that case, the electron orbital motion in a 
plane perpendicular to the axis of the dot is characterized 
by the same quantum number both in the dot and in the 
leads |l^ , and two orbitals participate in the S — T tran- 
sitions. The problem can then be mapped onto a special 
version of the two-impurity Kondo model. 

Following I p^ we now apply the " poor man scaling ap- 
proach" |2^] to the Hamiltonian H (|l^). Neglecting the 
irrelevant potential scattering phase shift [|T] and using 
the above mentioned cutoff procedure, a system of scaling 
equations is obtained, (cf. p7[|) 

dji/dlnd = - [(ji)2 + (j2)2] , dj2/dlnd = -2hj2 (15) 

(here ji ~ poJ^, j2 = PoJ^'^, d — pqD). The correspond- 
ing RG flow diagram has the fixed point ji = oo, but the 
resulting Kondo temperature Tk{5) turns out to be a 
sharp function of 5 . It is maximal when the T, S 
states are quasi degenerate, 5 <C Tk{5). The scaling in 
this case is governed by the effective integral j+ = ji -\-j2, 
and the system (15) is reduced to a single equation 



d3+/d\nd^-{j+f 



(16) 



with Tko = I?exp(— l/j+). In the opposite limit 5 3> 
Tk{5) the scaling of J'^^ stops &t D ~5. Then ji,2('5) = 

ji,2ln-i (^f^) and Tk{5) = ^exp [-l/ji(^)] « Tko- 
The singlet ground state S with zero Tk is realized when 
~5 <0,\E\>Tk(5). 

4. The salient features of FM stem from the qualitative 
dependence of its ground state and low-energy spectrum 
on the coupling constants V and Wi . The unusual singlet- 
triplet crossing should show up in the magnetic proper- 
ties of adsorbed molecules and tunnel transparency of 
asymmetric double quantum dots. 

According to quantum chemical calculations of the en- 
ergy spectrum of isolated cerocene molecule, the Van 
Vleck paramagnetic contribution of S" — T excitations 
is too weak to overcome the Larmor diamagnetic con- 
tribution of CsHg rings ||7|,pl|. This situation can dras- 
tically change for a FM adsorbed on a metallic layer. 
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The fixed point ji = oo corresponds to the scattering 
phases rjcriep) = 7''/2. In the case of adsorbed FM this 
means that the molecule has a residual spin 1/2 which 
interacts ferromagnetically with the conduction electrons 
p6[ . The temperature dependence of magnetic suscepti- 
bility x(2^) is predetermined by the energy parameters S 
and Tk{S). In particular, x(^) conserves its Curie-like 
character down to the lowest temperatures when S < 0, 
\5\ » Tk Then at T < the underscreened FM re- 
mains paramagnetic, and its susceptibility is 



x{T)^Xo{T)[1-Z{T/Tk)] 



(17) 



Here xo — iC /AT, C = {g^b)^, and Z{x) is the invari- 
ant coupling function (solution of the Gell-Mann - Low 
equation, see |2^]). The triplet spin state is restored at 
T > Tk . In this regime the Kondo corrections as well 
as admixture of singlet state can be calculated by per- 
turbation theory, with the result. 



X(T) 
xiT) 



2C[3 - exp(-Vr)] 
3T 
2C 

T[3 + cxp{-S/T)] 



1 - 1/ln 



T , D 



1 - 1/ln 



T 

Tko 



(18) 



respectively for S :s> Tko and S <C Tko- In the case of 
artificial FM the resonance scattering phase means per- 
fect tunneling transparency of the quantum dot at T = 
and a logarithmic fall off at high temperatures. To cal- 
culate the tunneling transparency of FM sandwiched be- 
tween two leads, one should add an index n = L, R to 
the operator Cnka- and switch to the standing wave basis 
V2cka± = CLka ± CRka (in & Symmetric configura- 
tion WiL — Win). Then only the wave (-I-) is involved 
in tunneling, and the zero bias anomaly in the differen- 
tial conductance G(T) (due to Kondo cotunneling) in the 
weak coupling regime T > Tk is found as in ([l8|): 



G/Go = 2ln-^[T/TK{6)]+.njiHD/S), 
G/Go = 3lii-^[T/TKo] 



(19) 



respectively for the two limiting cases (5 3> T 3> Tks and 
T ^ Tko ^ Here Go = Ane'^/h. Again the maximum 
effect is achieved in a nearly degenerate case. At T — > 
the conductance tends to the unitarity limit. 

In conclusion, the interplay between ZR-type coupling 
in real and/or artificial molecules and Kondo coupling 
between molecules and metallic reservoir may result in 
a crossover from a singlet spin state in a weak-coupling 
regime to an underscreened 5 = 1 state at zero T. The 
onset of Kondo regime in double quantum dot with even 
occupation can be driven either by a magnetic field or by 
a gate voltage. 
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